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Abstract 

We prove the most general theorem about spectral stability of multi-site breathers in the 
discrete Klein-Gordon equation with a small coupling constant. In the anti-continuum limit, 
multi-site breathers represent excited oscillations at different sites of the lattice separated by 
a number of "holes" (sites at rest). The theorem describes how the stability or instability 
of a multi-site breather depends on the phase difference and distance between the excited 
oscillators. Previously, only multi-site breathers with adjacent excited sites were considered 
within the first-order perturbation theory. We show that the stability of multi-site breathers 
with one-site holes change for large-amplitude oscillations in soft nonlinear potentials. We 
also discover and study a symmetry-breaking (pitchfork) bifurcation of one-site and multi-site 
breathers in soft quartic potentials near the points of 1:3 resonance. 

1 Introduction 

Space-localized and time-periodic breathers in nonlinear Hamiltonian lattices have been studied 
extensively in the past twenty years. These breathers model the oscillatory dynamics of particles 
due to the external forces and the interaction with other particles. A particularly simple model 
is the Klein-Gordon lattice, which is expressed by the discrete Klein-Gordon equation, 

iin + V'{Un) = €{Un+l - 2li„ + Un-l), n £ Z, (1) 

where t E M is the evolution time. Unit) G M is the displacement of the n-th particle, 1/ : M — )• M 
is a smooth on-site potential for the external forces, and e G M is the coupling constant of the 
linear inter-particle interaction. For the sake of clarity, we will assume that the potential V is 
symmetric, but a generalization can be formulated for non-symmetric potentials V. We will also 
assume that V'{u) can be expanded in the power series near w = by 

V'{u) = u±u^ + 0{u^) as n^O. (2) 

The plus and minus signs are referred to as the hard and soft potentials, respectively. 

A simplification of analysis of the Klein-Gordon lattice was proposed by MacKay and Aubry [7] 
in the anti-continuum limit of small coupling constant e — )• 0. This limit inspired many researchers 
to study existence, stability, and global dynamics of space-localized and time-periodic breathers 
[2]. Since all oscillators are uncoupled at e = 0, one can construct time-periodic breathers localized 
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at different sites of the lattice (called the multi-site breathers) . All these time-periodic multi-site 
breathers are uniquely continued with respect to the (small) coupling constant e if the period of 
these oscillations is identical and the oscillations are synchronized either in phase or out of phase. 

Spectral stability of discrete breathers, which are continued from the anti-continuum limit 
e = 0, was considered by Morgante et al. [8] with the help of numerical computations. These nu- 
merical computations suggested that spectral stability of small-amplitude breathers in the discrete 
Klein-Gordon equation ([1]) is similar to the spectral stability of solitons in the discrete nonlinear 
Schrdinger (DNLS) equation. 

The DNLS approximation for small-amplitude and slowly varying oscillations relies on the 
asymptotic solution, 

Unit) = [a„(et)e^* + a„(et)e-'*] + Oi^ (e^/^), (3) 

where e > is assumed to be small, r = et is the slow time, and {a„(T)}„gz E are amplitudes 
of harmonic oscillations with the linear frequency u = 1. Substitution of ([3l) to ([1]) yields the 
DNLS equation to the leading order in e, 

2id„ = On+i - 2a„ an_i =F 3|a„|^a„, n G Z. (4) 

The hard and soft potentials ([2]) result in the defocusing and focusing cubic nonlinearities of the 
DNLS equation dH), respectively. Existence and continuous approximations of small-amplitude 
breathers in the discrete Klein-Gordon equation were considered recently by Bambusi et al. [3]. 

Numerical results of [8] can be summarized as follows. In the case of the focusing nonlinearity, 
the only stable discrete soliton of the DNLS equation corresponds to the out-of-phase oscilla- 
tions of the excited sites of the lattice. This conclusion does not depend on the number of "holes " 
(oscillators at rest) between the excited sites at the anti-continuum limit. The stable oscillations 
in the case of the defocusing nonlinearity can be recovered from the stable out-of-phase oscillations 
supported by the focusing nonlinearity using the staggering transformation, 

a„(T) = (-l)"6„(T)e2-, (5) 

which changes the DNLS equation (|4]) to the form, 

2ibn = bn+i-2bn-\-bn-i±3\bn\'^bn, n £ Z. (6) 

Consequently, we have the following statement. In the case of the defocusing nonlinearity, the 
only stable discrete soliton of the DNLS equation with adjacent excited sites corresponds to the 
in-phase oscillations of the excited sites. 

The numerical observations of [8] were rigorously proved for the DNLS equation (jH by Peli- 
novsky et al. [9\ . Similar conclusions on the spectral stability of breathers in the discrete Klein- 
Gordon equation ([T|) were also reported in the literature under some simplifying assumptions. 

Archilla et al. [1\ used the perturbation theory for spectral bands to consider two-site, three- 
site, and generally multi-site breathers in Klein-Gordon lattices. Theorem 6 in [1] states that 
in-phase multi-site breathers are stable for hard potentials and out-of-phase breathers are stable 
for soft potentials for e > 0. The statement of this result misses however that the corresponding 



2 



computations are justified for multi-site breathers with adjacent excited sites: no "holes" (oscilla- 
tors at rest) in the limiting configuration at e = are allowed. More recently, Koukouloyannis and 
Kevrekidis [6] recovered exactly the same conclusion using the averaging theory for Hamiltonian 
systems in action-angle variables. To justify the use of the first-order perturbation theory, the 
multi-site breathers were considered to have adjacent excited sites and no holes. The equivalence 
between these two approaches was addressed by Cuevas et al. [5]. 

It is the goal of our paper to rigorously prove an analogue of this statement for all multi- 
site breathers, including breathers with holes between excited oscillators in the anti-continuum 
limit. We will use perturbative arguments for characteristic exponents of the Floquet monodromy 
matrices. To be able to work with the higher-order perturbation theory, we will combine these 
perturbative arguments with the theory of tail-to-tail interactions of individual breathers in lattice 
differential equations. Although the tail-to-tail interaction theory is well-known for continuous 
partial differential equations [TI] , it is the first time to our knowledge when this theory is extended 
to nonlinear lattices. 

Multi-site breathers with holes have been recently considered by Yoshimura [.15j in the context 
of the diatomic Fermi-Pasta-Ulam lattice near the anti-continuum limit. In order to separate 
variables n and t and to perform computations using the discrete Sturm theorem (similar to the 
one used in the context of NLS lattices in |9j), the interaction potential was assumed to be nearly 
homogeneous of degree four and higher. Similar work was performed for the Klein-Gordon lattices 
with a purely anharmonic interaction potential [16]. Compared to this work, our treatment is 
valid for a non-homogeneous on-site potential V satisfying expansion Q and for the quadratic 
interaction potential. 

Within our work, we have discovered new important details on the spectral stability of multi- 
site breathers, which were missed in the previous works [H EJ [8]. In the case of soft potentials, 
breathers of the discrete Klein-Gordon equation ([1]) can not be continued far away from the 
small-amplitude limit described by the DNLS equation (|4} because of the resonances between 
the nonlinear oscillators at the excited sites and the linear oscillators at the sites at rest. These 
resonances change the stability conclusion. In particular, the out-of-phase oscillations may become 
unstable in soft nonlinear potentials even if the coupling constant is sufficiently small. 

Another interesting feature of soft potentials is the symmetry-breaking (pitchfork) bifurcation 
of one-site and multi-site breathers that occur near the point of resonances. In symmetric po- 
tentials, the first non-trivial resonance occurs near w = ^, that is, at 1:3 resonance. We analyze 
this bifurcation by using asymptotic expansions and normal form equations. It is interesting that 
the normal form equation for 1:3 resonance which we analyze is different from the normal form 
equations considered in the previous studies of 1:3 resonance [H [T2| I13j . because we are looking 
at bifurcations of periodic solutions far from the equilibrium points. 

The paper is organized as follows. Existence of space-localized and time-periodic breathers at 
the anti-continuum limit is reviewed in Section 2 using a new version of the tail-to-tail interaction 
theory. The main result on spectral stability of multi-site breathers for small coupling constants 
is formulated and proved in Section 3. Section 4 compares Floquet theory with the Hamiltonian 
averaging theory to show equivalence of our conclusions with those reported earlier in [H [6]. 
Section 5 illustrates the spectral stability of breathers in soft potentials near the resonances 
numerically. Section 6 is devoted to studies of the symmetry-breaking (pitchfork) bifurcation 
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using asymptotic expansions and normal forms for the 1:3 resonance. Section 7 summarizes our 
findings. 



2 Existence of multi-site breathers near the anti-continuum limit 

In what follows, we will use bold-faced notations for vectors in discrete space /*(Z) defined by 
their norms 

\\u\\is := I ^ Im^I'' J , s > 1. 

Components of u are denoted by Un for n £ Z. These components can be functions of t, in which 
case they are considered in Hilbert-Sobolev spaces Hp^^(0,T) of T-periodic functions for s > 0, 
equipped with the norm, 
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where {cmjmez is the set of Fourier coefficients in the Fourier series of a T-periodic function /, 

N f 2TTimt\ 

/(*)= 5^c^expf , t€[0,T]. 

We consider space-localized and time-periodic breathers u G Hp^j.(0,T)) of the discrete 
Klein-Gordon equation ([T|) with smooth even V and e > 0. Parameter T > represents the 
fundamental period of the time-periodic breathers. Accounting for symmetries, we shall work in 
the restriction of Hp^j.{0,T) to the space of even T-periodic functions, 

Hiio,T) = {f e h;^MT) : f{-t) = f{t), tern}, s>o. 

At e = 0, we have an arbitrary family of multi-site breathers, 

uW(t) = j;afc^(t)e,, (7) 

kes 

where is the unit vector in /^(Z), S" C Z is a finite set of excited sites of the lattice, G {+1, — 1} 
encodes the phase factor of the k-th. oscillator, and f G i^pgj,(0,T) is a solution of the nonlinear 
oscillator equation at the energy level E, 

^ + V'i^) = E = ^^^ + Vi^). (8) 

Period of oscillations T is uniquely defined by the energy level E, 

T = V2r , , (9) 

7-a ^E - V[ip) 

where a is the smallest positive root of V{a) = E. 
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Figure 1: Top: the phase plane {^p,^) for the hard (left) and soft (right) potentials (llOp . Bottom: 
the period T versus energy E for the hard (solid) and soft (dashed) potentials. 



Remark 1 All nonlinear oscillators at the excited sites of S C Z in the limiting configuration 
^ have the same period T. Two oscillators at the j-th and k-th sites are said to he in phase if 
GjCTj. = 1 and out of phase if CTjaj. = —1. 

Our main example of the nonlinear potential V is the truncated expansion ([2|) with 

V'{u) = u±u\ (10) 

The dependence T versus E is computed numerically from ([9]) and is shown on Figure [2] together 
with the phase portraits of the system ([8]). For the hard potential with the plus sign, the period 
T is a decreasing function of E in (0, 27r), whereas for the soft potential with the minus sign, the 
period T is an increasing function of E with T > In. 



Persistence of the limiting configuration ([7]) as a space-localized and time-periodic breather of 
the discrete Klein-Gordon equation ([T|) for small values of e is established by MacKay & Aubry 
[7]. The following theorem gives the relevant details of the theory that are useful in our analysis. 
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Theorem 1 Fix the period T and the solution ip G H^{0,T) of the nonlinear oscillator equation 
^ with an even V G C°°(M) satisfying ^ and assume that T 7^ 2'Kn, n S N and T'{E) ^ 0. 
Define u^^^ by the representation ([^j with fixed S d'L and {(Tk}k&s- There are cq > and C > 
such that for all e G (— eoj^o); there exists a unique solution vS-''^ G P{7j, H'^{0,T)) of the discrete 
Klein-Gordon equation (OP satisfying 

||u(') - u(°)||;2(2_^2^^(o,T)) < C\e\. 

Moreover, the map R3 u(') G P{Z, H'^{0, T)) is C°° at e = 0. 

Proof. Thanks to the translational invariance of the nonhnear oscillator equation ([8]), we can 
fix an even ip and consider H^(0,T), the restriction of Hp^^{0,T) to even functions. Under the 
condition T'{E) 7^ 0, operator 

L, = 92 + V'i^it)) : i/2(o,T) ^ LliO,T) 

is invertible, because the only eigenvector ip oi L = + V"{ip{t)) : H^^^{Q,T) — )• Lpgj,(0, T) is 
odd in t. Similarly, operator 

Lo = 92 + 1 : i/2,,(0, T) ^ L2g,(0, T) 

is invertible if T 7^ 2im, n G N. 

Substituting u = u*^^^ +w, where u(0) is the limiting breather ([7]) to the discrete Klein-Gordon 
equation ([1]), we obtain the coupled system of differential-difference equations 

LeWn = erfji - 2u(f) + uf_^) + Af„(w, e), n G 5 

and 

1+1 



Lom„ = e{u^^l^ - 2u^^^ + u^^l^) + Ar„(w, e), n G Z\S, 



where 

Thanks to Banach algebra of and the smoothness of V , we have N(vir, e) : xR^ H^. For 
any w in a ball in i/g centered at with radius J > 0, there are constants Cs,Ds > such that 



and 

||N(wi,e) - N(v^^2,e)||/2(2,H2,^(o,T)) < -D^ (^e + ||wi ||i2(a_/^2^^(o,T)) + l|w2||p(z,/f2^^.(o,T)) 

X ||vi^i - W2||/2(2_^2^^,(o,T))- 

Thanks to the invertibility of the linearized operators Lg and Lq on Lg, the result of the theorem 
follows from the Implicit Function Theorem. Because V G C°°(M), the map e i— >■ u^*^) is C°° at 
e = 0. I 
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Remark 2 Although persistence of other breather configurations, where oscillators are neither in 
phase nor out of phase, can not be apriori excluded, we restrict our studies to the most important 
and physically relevant breather configurations covered by Theorem [7J 

We shall now introduce the concept of the fundamental breather for the set S = {0}. Multi- 
site breathers for small e > can be approximated by the superposition of fundamental breathers 
at a generic set S of excited sites up to and including the order, where the tail-to-tail interactions 
of these breathers occur. 

Definition 1 Letu^^^ G H'^{0,T)) be the solution of the discrete Klein-Gordon equation 
for small e > defined by TheoremlJl for a given u^^\t) = (p{t)eQ. This solution is called the 
fundamental breather and we denote it by cf)^^\ 

By Theorem [H we can use the Taylor series expansion, 

A:=0 ■ ^=0 

up to any integer A'^ > 0. Thanks to the discrete translational invariance of the lattice, the 
fundamental breather can be centered at any site j G Z. Let tj : /^(Z) — t- /^(Z) be the shift 
operator defined by 

{Tju)n = Un-j, n G Z. 

If 0*-"^^ is centered at site 0, then Tjcf)^^^ is centered at site j G Z. The simplest multi-site breather 
is given by the two excited nodes at j G Z and k £ 1^ with j ^ k. 

Lemma 1 Let u^^\t) = aj(p{t)ej + (7i^Lp{t)ek with j ^ k and = |j — A;| > 1. Let u^^^ G 
l'^{Z, H^[0,T)) be the corresponding solution of the discrete Klein-Gordon equation (CP for small 
e > defined by TheoremUl Let {ipm}m=i ^ ^e(0'^) defined recursively by 

Loipm := {dt +l)ipm = V>m-i, m = l,2,...,N, (12) 
starting with (po = ip, and let ip]\f G H^{0,T) be defined by 

Le^N := (d^ + V'i^m^PN = (13) 

Then, we have 

u(^) = a,-r,0(^'^) + afcTfc0(^'^) + (a.-e^ + a^e,-) (^^v - ^n) + ©/^(^^^^^^^(o.t)) (e^+'). (14) 



Proof. By Theorem[Tl the limiting configuration u('^)(t) = aj(p{t)ej +ai^Lp{t)ei^ with two excited 
sites generates a C°° map, which can be expanded up to the + 1-order, 



u 



(^) 



N 

E 



d'' 

'k\d? 



u 



(^) 



e=0 



+ Cp(Z,J/2,,(0,T))(e^^^; 



(15) 
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Substituting (|15p into ([T]) generates a sequence of equations at each order of e, which we consider 
up to and including the terms of order N. 

The central excited site at n = in the fundamental breather (f)^''^ generates fluxes, which reach 
sites n = ibm at the m-th order. Because t/)^*^'™^ is compactly supported on {—m, —m + 1, m} 
and all sites with n 7^ contain no oscillations at the 0-th order, we have 

= e™^^, (16) 

where {ipm}m=i ^ -f^e(O)^) computed from the linear inhomogeneous equations (|T2]l starting 
with (fQ = (p. Note that equations (|12p are uniquely solvable because T 7^ 2Tm, n G N. 

For deflniteness, let us assume that j = and k = N > 1. The fluxes from the excited sites 
n = and n = N meet at the A^/2-th order at the middle site n = N/2 if A'^ is even or they overlap 
at the {N + 1) /2-th order at the two sites n = (iV - l)/2 and n= {N + l)/2\i N is odd. In either 
case, because of the expansion ([2|), the nonlinear superposition of these fluxes affects terms at 
the order 3A^/2-th or 3(A'^ + l)/2-th orders, that is, beyond the A^-th order of the expansion (fT4l) . 
Therefore, the nonlinear superposition of fluxes in higher orders of e will deflnitely be beyond the 
A^-th order of the expansion (|14p . 

Up to the A^-th order, all correction terms are combined together as a sum of correction terms 
from the decomposition (jlip centered at the j-th and A;-th sites, that is, we have 

u(e) = cT,r,<^(^'^-i)(e) + afcrfc0(^'^-i)(e) + Op(z,//2^,(o,T))(e^). (17) 

At the A'^-th order, the flux from j-th site arrives to the fc-th site and vice versa. Therefore, 
besides the A'^-th order correction terms from the decomposition (jlip . we have additional terms 
(cTj-efc + akGj) i/j^ at the sites n = j and n = k. Thanks to the linear superposition principle, 
these additional terms are given by solutions of the inhomogeneous equations p3|) . which are 
uniquely solvable in H^{0,T) because T'{E) ^ 0. We also have to subtract {o'jGk + ^k^j) 
from the A'^-th order of ajTjip^^'^^ + akTk(i)^'''^\ because these terms were computed under the 
assumption that the /c-th site contained no oscillations at the order for (TjTjcl)^''^^ and vice versa. 
Combined all together, the expansion is justified up to terms of the A'-th order. ■ 

3 Stability of multi-site breathers near the anti-continuum limit 

Let u G H^{0,T)) be a multi-site breather in Theorem [1] and e > be a small parameter 
of the discrete Klein-Gordon equation ([T]). When we study stability of breathers, we understand 
the spectral stability, which is associated with the linearization of the discrete Klein-Gordon 
equation ([TJ by using a perturbation w(t) in the decomposition u(t) -|-w(t). Neglecting quadratic 
and higher-order terms in w, we obtain the linearized discrete Klein-Gordon equation, 

Wn + V"{Un)Wn = e{Wn+l - 2Wn + tt^n-l), n € Z. (18) 

Because u(t + r) = u(t), an infinite-dimensional analogue of the Floquet theorem applies and the 
Floquet monodromy matrix Ai is defined by w(T) = 7Ww(0). We say that the breather is stable 
if all eigenvalues of A^, called Floquet multipliers, are located on the unit circle and it is unstable 
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if there is at least one Floquet multiplier outside the unit disk. By reversibility of the linearized 
system (llSp in t, Floquet multipliers come in pairs m and ^2 with fiiiJ,2 = 1- 

For e = 0, Floquet multipliers can be computed explicitly because A4 is decoupled into a 
diagonal combination of 2-by-2 matrices {Mn}nez, which are computed from solutions of the 
linearized equations 

w^ + Wn = 0, n G Z\S (19) 

and 

iVn + V"{if)Wn = 0, n G S. 
The first problem (I19p admits the exact solution, 



(20) 



Wn{t) = Wn{0) COs(t) + Wn{0) sin(t) 



Mr, 



cos(r) sin(r) 
— sin(r) cos(T) 



n e Z\S. (21) 

Each Mn for n G Z,\S has two Floquet multipliers at fii = e*^ and fi2 = e~*"^- If T 7^ lirn, n G Z, 
the Floquet multipliers are located on the unit circle bounded away from the unit multiplier. 
The second problem (I20p also admits the exact solution, 

».W = ^^W + ^8.,(*). n.S, (22) 

where (p(t) is a solution of the nonlinear oscillator equation ([8]) with the initial condition 

(/j(0) = a, ip{0) = 0. 

Since a is the smallest positive root of V{a) = E for a fixed value of E and ip{t) is T-periodic, we 
have 



(p{Q) = -V'{a), dE^{T) = a'{E) 
As a result, we obtain, 



1 



V'{a) 



dE^{T) = -(p{T)T'{E) = V'{a)T'{E). 



Mr, 



n£S. 



1 

r{E){V'{a)f 1 

Note that V'{a) / (or T is infinite). If T'{E) / 0, each M„ for n G 5 has the unit Floquet 
multiplier of geometric multiplicity one and algebraic multiplicity two. 

We conclude that if T 7^ 2-iTn, n G N and T'{E) 7^ 0, the limiting multi-site breather ([7]) at the 
anti-continuum limit e = has an infinite number of semi-simple Floquet multipliers at /^i = e*^ 
and //2 = e~*"^ bounded away from the unit Floquet multiplier of algebraic multiplicity 2|S'| and 
geometric multiplicity \S\. 

Semi-simple multipliers on the unit circle are structurally stable in reversible dynamical sys- 
tems. Therefore, the instability of the multi-site breather may only arise from the splitting of the 
multiple unit multiplier for e 7^ 0. 

To consider Floquet multipliers, we can introduce the characteristic exponent A in the decom- 
position w(t) = W(t)e-^*. If /u = e^^ is the Floquet multiplier of the monodromy operator Ai, 
then W G /^(Z, //pgj,(0, T)) is a solution of the eigenvalue problem. 



Wr, + V"iUn)Wn + 2XWn + X^Wn = e{Wn+l " 2Wn + IVn-l), n G 



(23) 
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In particular, Floquet multiplier ^ = 1 corresponds to the characteristic exponent A = 0. The 
generalized eigenvector Z of the eigenvalue problem (j23p for A = solves the inhomogeneous 
problem, 

Zn + V"(Un)Zn = e{Zn+l-2Zn + Zn-l)-2Wn, n G Z, (24) 

where W is the eigenvector of (I23p for A = 0. 

At e = 0, the eigenvector W of the eigenvalue problem (I23p for A = is spanned by the linear 
combination of |5| fundamental solutions, 

MV^^){t) = Y,cw{t)eu. (25) 

k&S 

The generalized eigenvector Z is spanned by the linear combination of \S\ fundamental solutions, 

Z(°Hi) = -Y. Cfe^Wefc, V := 21'^ Cp, (26) 
fees' 

where = + V"{tp{t)) : H'^{0,T) LI{0,T) is invertible and ^3 G LI{0,T) (see the proof of 
Theorem [T|) . 

Because of the translational invariance in t, we note that if u = cf)^''^ is the fundamental 
breather in Definition [U then W = dt4)^^^ = O^""^ is the eigenvector of the eigenvalue problem (j23p 
for A = and small e > and there exists a generalized eigenvector Z = fi^^^ of the inhomogeneous 
problem (|24p . which exists because dtO^^^ has the opposite parity in t compared to 6^''\ 

By Taylor series expansions (llip . for any integer > 0, we have 



where 0^^'^^ and fi^'''^'> are polynomials in e of degree N. It follows from equations (|25p and (j26p 
that 

0W = Vi(Oeo, /xW = -t;(t)eo. (28) 

This formalism sets up the scene for the perturbation theory, which is used to prove the main 
result on spectral stability of multi-site breathers. We start with a simple multi-site breather 
configuration with equal distances between excited sites and then upgrade this result to multi-site 
breathers with non-equal distances between excited sites. 

Lemma 2 Let u(°)(t) = T.jLiCrjtf{t)ejN with fixed M,N G N. Let u(^) G P{Z,H^{0,T)) be 
the corresponding solution of the discrete Klein-Gordon equation (OP for small e > defined by 
TheoremUl Let {ipm}m=o defined by Lemma ll\ starting with cpo = if. Then the eigenvalue 
problem i23\) for small e > has 2M small eigenvalues, 



A = e^/2A + 0(e(^+i)/2), 
where A is an eigenvalue of the reduced eigenvalue problem 



^P^K^c = KnMc, cgC^. (29) 
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Here ^ 

Kn = (f(fN-idt 
Jo 

and the matrix 7W G M^^^*^ is defined by its entries, 

Mij = -aj{aj-i + aj+i)6ij + 6ij+i + dij-i, i,j G {1,2, ...,M}, ctq = ctaz+i = 0, 
where 6ij is the Kronecker symbol. 

Proof. At e = 0, the eigenvalue problem (j23p admits eigenvalue A = of geometric multiplicity 
M and algebraic multiplicity 2M, which is isolated from the rest of the spectrum. Perturbation 
theory in e applies thanks to the smoothness of u*^*^-* in e and V' in u. Perturbation expansions 
(so-called Puiseux series, see recent work [14J) are smooth in powers of e^/^ thanks to the Jordan 
block decomposition at e = 0. Therefore, we are looking for the eigenvectors of the eigenvalue 
problem (|23p in the subspace associated with the eigenvalue A = using the decompositions, 

A = e^/^X (30) 

and 

M 

i=i 

M 

+6^/2AX^c,W^'^*) + e^W, (31) 
i=i 

where N^, = N/2 if is even and = {N — l)/2 if N is odd, whereas W is the remainder term 
at the A^-th order in e. The decomposition formula (I3ip follows from the superposition (I14p up to 



the A^-th order in e. The fluxes from '^ji^{j-i)N + ^{j+i)N)''pN in the superposition ([T 

are to be accounted at the equation for W. Note that our convention in writing (|3ip is to drop 
the boundary terms with Goat and e(^]^,j^i^]\f. 

Substituting (|3Up and (|3ip to (|23p . all equations are satisfied up to the A^-th order. At the A^-th 
order, we divide (|23p by and collect equations at the excited sites n = jN for j G {1,2,..., M}, 

WjN + V"{ip)WjN = (cj+i + Cj-i)ipN-i - crj{(^j+i + crj-i)cjV"'{ip)'ipN^ 

+ ~X^Cj{2v-ip) + 0{e^/^), (32) 

where we admit another convention that (Tq = (Tm+i = and cq = cm+i = 0. To justify (j32p . the 
term (pN-i comes from the fluxes from n = {j + 1)N and n = [j — 1)N sites generated by the 
derivatives of the linear inhomogeneous equations (|12p and the term aj{aj+i+aj-i)cjV"'{(p)'ipNip 
comes from the expansion (|14p of the nonlinear potential V"{ujN). 

Expanding A = A + 0{e^^'^) and projecting the system of linear inhomogeneous equations (j32p 



to 99 G H^^^{0,T), the kernel of L = d'^ + V"{ip) : Hl^^{0,T) Ll^^{0,T), we obtain the system 
of difference equations, 

A^Cj / [ip"^ + 2v(p) dt = {cj+i + Cj-i) / ifipN-idt - aj{aj^i + aj-i)cj / V'" {ip)ipNV''^dt, 
Jo Jo Jo 
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where the integration by parts is used to simphfy the left-hand side. Differentiating the hnear 
inhomogeneous equation (fT3|l and projecting it to ip, we infer that 

Jo Jo 
The lemma is proved if we can verify that 

To do so, we note that v = 2L~^Cp is even in t G M, so that it is given by the exact solution, 

where C G M. From the condition of T-periodicity for v{t) and i'(t), we obtain 

y(O) = v(T) = Ca'{E), v{0) = = v{T) = Tip{0) - CT'{E)if{0), 
hence C = T{E)/T'{E) and 

/ {ip^ + 2v(p) dt = 2C (pdE^dt = -C {^BeP + V'{ip)dEv) dt 
Jo Jo Jo 

where equation ([8|) has been used. The proof of the lemma is complete. ■ 

We shall now count eigenvalues of the matrix eigenvalue problem (I29p to classify stable and 
unstable configurations of multi-site breathers near the anti-continuum limit. 

Lemma 3 LetuQ be the numbers of negative elements in the sequence {(Tj(Tj+i}^Y^- IfT'{E)Ki^ > 
0, the matrix eigenvalue problem \2y\) has exactly uq pairs of purely imaginary eigenvalues A and 
M — I — uq pairs of purely real eigenvalues fi counting their multiplicities, in addition to the double 
zero eigenvalue. If T' {E)Kn < 0, the conclusion changes to the opposite. 

Proof. We shall prove that the matrix M has exactly no positive and M — 1 — uq negative 
eigenvalues counting their multiplicities, in addition to the simple zero eigenvalue. If this is the 
case, the assertion of the lemma follows from the correspondence A-^ = — '^j,f|^^^ 7, where 7 is an 
eigenvalue of M. 

Setting Cj = CTjbj, we rewrite the eigenvalue problem Mc = 7c as the difference equation, 

- bj) + ajaj-iibj-i - bj) = 76^-, j G {1, 2, M}, (33) 

subject to the conditions ao = o"m+i = 0. Therefore, 7 = is always an eigenvalue with the 
eigenvector b = 1 = (1,1, The coefficient matrix in (|33p coincides with the one analyzed 

by Sandstede in Lemma 5.4 and Appendix C |10j. This correspondence yields the assertion on 
the number of eigenvalues of A^. ■ 
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Before formalizing the results of Lemmas [2] and [3] to general multi-site breathers, we consider 
two examples, which are related to the truncated potential (llOp . Let us use the Fourier cosine 
series for the solution ip G H^{0,T), 

<^W = 5^c„(i?)cos('^V (34) 

for some square summable set {c„(£')}„gN- Because of the symmetry of V, we have ip{T/4) = 0, 
which imply that Cn{E) = for all even n G N. Solving the linear inhomogeneous equations (|12p . 
we obtain 

Using Parseval's equality, we compute the constant Kj^ in Lemma [2l 

r . ^ , . 2 V- T2^-=^(^)n2|c„(^)|2 



(36) 



For the hard potential with V'{u) = u + v? , we know from Figure [2] that T{E) is a decreasing 
function from r(0) = 2ti to \miE^ocT{E) = 0. Since T'{E) < and T{E) < 27r, we conclude 
that T'{E)Kn < if TV is odd and T'{E)Kn > if TV is even. By Lemma El if N is odd, the 
only stable configuration of the multi-site breathers is the one with all equal {aj}j^i (in-phase 
breathers), whereas if is even, the only stable configuration of the multi-site breathers is the 
one with all alternating {cjj}*£^ (out-of-phase breathers). This conclusion is shown in the first 
fine of Table L 

For the soft potential with V'{u) = u — u^, we know from Figure [2] that T{E) is an increasing 
function from T(0) = 27r to lim^;^^;;, T{E) = oo, where Eq = ^. If T{E) is close to but above 27r, 
then the first term in the series (j36p dominates and > for all N G N. At the same time, 
T'{E) > and Lemma [3] implies that the only stable configuration of the multi-site breathers 
is the one with all alternating {(Tj}^^ (out-of-phase breathers). The conclusion holds for any 
T > 27r if is odd, because > in this case. 

This precise conclusion is obtained in Theorem 3.6 of |9] in the framework of the DNLS 
equation It is also in agreement with perturbative arguments in [1] [6j, which are valid for 
A^ = 1 (all excited sites are adjacent on the lattice). To elaborate this point further, we show 
in Section H] the equivalence between the reduced eigenvalue problem (I29p with A^ = 1 and the 
criteria used in [1] |6]. 

For even A^ S N, we observe a new phenomenon, which arise for the soft potentials with larger 
values of T[E) > 27r. The constant K]\f in ([36|) may change sign to < 0, in which case, the 
conclusion on stability of the multi-site breather change directly to the opposite: the only stable 
configuration of the multi-site breathers is the one with all equal {fTj}^'£^ (in-phase breathers). In 
particular, this happens when T{E) is close but below dn, which corresponds to the first resonance 
1:3. This resonance is further studied in Sections [5] and |6j This conclusion is shown in the second 
line of Table I. 
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N odd 


N even 


hard potential 
V'{u) = u + u^ 


in-phase 


anti-phase 


soft potential 
V'{u) = u — u^ 


anti-phase 


2Tr < T < anti-phase 
< T < Gvr in-phase 



Table I: Stable multi-site breathers in hard and soft potentials. The stability threshold 
corresponds to the zero of Ki\f. 

We conclude this section with the stability theorem for general multi-site breathers. For the 
sake of clarity, we formulate the theorem in the case when T'(E) > and all Kjsf > 0, which 
arises for the soft potential with odd N. Using Lemma [3l the count can be adjusted to the cases 
of T'{E) < and/or Kn < 0. 

Theorem 2 Let {nj}jL^ C Z be an increasing sequence withM e N. Leiu(^) e P{Z,H^{0,T)) be 

a solution of the discrete Klein-Gordon equation (OP in Theorem\^with \i^^\t) = o'j'/?(t)e„^. 
for small e > 0. 

Define {A'j}^^^"'^ and {Kfq.}^ls[^ by Nj = Wj+i — nj and K]\f_. = ipip]\i._idt. Assume 
T'{E) > and Kn^ > for all Nj. 

Let uq be the numbers of negative elements in the sequence {crjCrj+i}*^^^- eigenvalue 
problem i23\) at the discrete breather u^^^ has exactly hq pairs of purely imaginary eigenvalues X 
and M — 1 — uq pairs of purely real eigenvalues A counting their multiplicities, in addition to the 
double zero eigenvalue. 

Proof. The limiting configuration u('^)(i) = Cj'/3(t)e„^. identifies clusters of excited sites 

with equal distances Nj between the adjacent excited sites. 

According to Lemma [21 splitting of M double Jordan blocks associated to the decompositions 
(|25p and (|26p occurs into different orders of the perturbation theory, which are determined by the 
set {Nj}jLl^. At each order of the perturbation theory, the splitting of eigenvalues associated with 
one cluster with equal distance between the adjacent excited sites obeys the matrix eigenvalue 
problem (I29p . which leaves exactly one double eigenvalue at zero and yields symmetric pairs of 
purely real or purely imaginary eigenvalues, in accordance to the count of Lemma [3j 

The double zero eigenvalue corresponds to the eigenvector W and the generalized eigenvector 
Z generated by the translational symmetry of the multi-site breather bifurcating from a partic- 
ular cluster of excited sites in the limiting configuration u^^\ The splitting of the double zero 
eigenvalues associated with the cluster happens at the higher orders in e, when the fiuxes from 
adjacent clusters reach each others. Since the end-point fiuxes from the multi-site breathers are 
equivalent to the fluxes (jl6p generated from the fundamental breathers, they still obey Lemma [1] 
and the splitting of double zero eigenvalue associated with different clusters still obeys Lemma [2j 

At the same time, the small pairs of real and imaginary eigenvalues arising at a particular order 
in e remain at the real and imaginary axes in higher orders of the perturbation theory because 
their geometric and algebraic multiplicity coincides, thanks to the fact that these eigenvalues are 
related to the eigenvalues of the symmetric matrix Ai in Lemma [3j 

Avoiding lengthy algebraic proofs, these arguments yield the assertion of the theorem. ■ 
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4 Comparison of Floquet theory with Hamiltonian averaging 



We will show here that the criterion of Lemma [3] for = 1 agrees exactly with the main conclu- 
sions of the previous studies [HIS] (see also [5j). 

The reduced eigenvalue problem (j29p for = 1 can be written in the form 



where 



A{E)K^c = Mc, c e C^^ 
AiE) - 



(37) 



(38) 



We will show that the quantity A(E) arises both in the spectral band theory used in [Tj and in 
the Hamiltonian averaging theory used in [6]. 

For the spectral band theory, we consider solutions of the spectral problem 



Lu = An, L = + V"{ip{t)) : HLMT) ^ LLMT). 



(39) 



Let M = $(r) be the monodromy matrix computed from the fundamental matrix solution <I>(t) 
of the system 



d 



dt 



(40) 



1 

X-V"{^{t)) 

subject to the initial condition ^*(0) = I £ M^^^. Since det(M) = 1, the Floquet multipliers /xi 
and fi2 satisfy 

Hifi2 = h ^Il + ^l2=tT{M) = F{X). 
In particular, = /i2 = 1 if tr(M) = 2, which is true at A = thanks to the exact solution (l22l) . 



a'{E) m 
a'(E) (piQ) J 



^>(r) 



1 

T'{E){V'{a)f 1 



Hence, we have -F(O) = 2. We will show that A{E) in (j38p determines the sign of -F'(O). Denote 
elements of M = $(r) by Mi^ for 1 < z, j < 2. Since det(M) = 1 for all A, we obtain 



F'(0) =aA(Mn+M22)|A=o 
Let U{t) be a solution of 



M2idxM_ 



12|A=0 



T'{E){V'ia)fdxM_ 



12|A=0- 



U{t) + V"{^{t))U{t) = ^{t), 



(41) 



subject to the initial condition U{0) = U{0) = 0. Then, U{T) = dxMi2\x=o'p{0). Solving the 
second-order equation (j4ip . we obtain an explicit solution 



U{t) = ^{t) ( 1 
from which we find that 



<^{s)dE'f{s)ds^ + dE<^{t) j (p'^{s)ds, 
F'(0) = -T'{E) r ip'^it)dt - ^^^^^ 

^0 



AiE)- 
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If T'{E) < (for the hard potentials), we have F'(0) > 0, which imphes that the spectral 
band of the purely continuous spectrum of operator L in L^(]R) is located to the left of A = 0. If 
T'{E) > (for the soft potentials), we have -F'(O) < 0, which implies that the spectral band of 
L in L^(R) is located to the right of A = 0. If A = uj{k) is the dispersion relation of the spectral 



band for k £ 



T(E) ' T{E) 



then near A = 0, we have 



uj{k) = -'^^^k'^ + 0{k^) as k^O uj"{0) = -2A{E). 

The identity a;"(0) = —2A{E) establishes the equivalence of the reduced eigenvalue problem (j37l 
with the spectral band theory considered in pLj. 

For the Hamiltonian averaging theory, we consider the action variable for the nonlinear oscil 
lator (HI), 

ME) , 

J = 4 / ^y2{E - V{ip))d^p. 
Jo 

Explicit computation shows that 

^ = 2V2 / , ^ = T(E). 

dE Jo ^E-V{ip) ' 

If cj = is the frequency of oscillations, then 

dE _ d^E _duj _ T'{E) _ 1 



dJ ' ' dp dJ T^E) T{E)A{E) ip^dt 

Therefore, the signs of ^(-E) and E"{J) coincide and this establishes the equivalence of the reduced 
eigenvalue problem (|37p with the Hamiltonian averaging theory considered in [6]. 



5 Numerical results 

We illustrate our analytical results on existence and stability of discrete breathers near the anti- 
continuum limit by using numerical approximations. The discrete Klein-Gordon equation ([T|) can 
be truncated at a finite system of differential equations by applying the Dirichlet conditions on a 
bounded subset of Z. 



5.1 Three-site model 

The simplest model which allows gaps in the initial configuration u^") is the one restricted to 
three lattice sites, e.g. n G { — 1,0, 1}. We choose the soft potential V'{u) = u — v? and rewrite 
the truncated discrete Klein-Gordon equation as a system of three Duffing oscillators with linear 
coupling terms, 

f uo + uo-ul = e{ui - 2uo + n_i), 
\ u±i + u±i-u^^ = e{uo-2u±i). 
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A fast and accurate approach to construct T-periodic solutions for this system is the shooting 
method. The idea is to find a E such that the solution u{t) E C^(]R+,M^) with initial 
conditions u(0) = a, u(0) = satisfy the conditions of T-periodicity, u(T) = a, u(r) = 0. 
However, these constraints would generate an over-determined system of equations on a. To set 
up the square system of equations, we can use the symmetry t i— )• — t of system (j42p . If we 
add the constraint u{T/2) = 0, then even solutions of system (I42p satisfy u(— T/2) = u{T/2) 
and u(— T/2) = —u{T/2) = 0. Periodicity of the resulting solution with period T then follows 
from translation invariance in t. Hence, the values of a E become the roots of the vector 
F(a) = u(T/2) E M^. 

We now construct a periodic solution u to system (j42p that corresponds to the anti-continuum 
limit u(0) as follows. Using the initial data u'-'^^(O) as an initial guess for the shooting method for 
a fixed value of T, we continue the initial displacement u(0) with respect to the coupling constant 
e > 0. After that, we fix a value of e and use the shooting method again to continue the initial 
displacement u(0) with respect to period T. 

Let us apply this method to determine initial conditions for the fundamental breather, 

4°) = </., nS = 0, (43) 

and for a two-site breather with a hole, 

nf = 0, u'il = ip. (44) 

In both cases, we can use the symmetry u_i(t) = ui{t) to reduce the dimension of the shooting 
method to two unknowns oq and ai. 

Figure [2] shows solution branches for these two breathers on the period-amplitude plane by 
plotting T versus oq and Oi. For 2tt <T < Gvr, solution branches are close to the limiting solutions 
(dotted line), in agreement with Theorem [TJ However, a new phenomenon emerges near T = Gvr: 
both breather solutions experience a pitchfork bifurcation and two more solution branches split off 
the main solution branch. The details of the pitchfork bifurcation occurring on the fundamental 
solution branch are shown on the insets of Figure [2j We checked numerically that the value of 
period T at the bifurcation point diverges to infinity in the limit of small e. 

Figure [2] also shows two branches of solutions for T > Gvr with negative values of oq and ai. 
One of the two branches is close to the breathers at the anti-continuum limit, as prescribed by 
Theorem [TJ We note that the breather solutions prescribed by Theorem [1] for T < Gvr and T > Gvr 
belong to different solution branches. 

Figure [3] shows the fundamental breather before (T = Svr) and after (T = S.Svr) pitchfork 
bifurcation. The symmetry condition u(T/4) = for the solution at the main branch is violated 
for two new solutions that bifurcate from the main branch. 

Let us now illustrate the stability result of Theorem [2] using the fundamental breather (j43p 
and the breather with a hole (I44p . We draw a conclusion on linearized stability of the breathers 
by testing whether their Floquet multipliers, found from the monodromy matrix associated with 
the linearization of system (I42p . stay on the unit circle. 
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Figure 2: The initial displacements ao and oi for the T-periodic solutions to system (I42p with 
e = 0.01. The solid and dashed lines correspond to the fundamental (|43p and two-site (j44l) 
breathers respectively. The dotted lines correspond to the T-periodic solutions to equation ([8]). 
The insets show the pitchfork bifurcation of the fundamental breather. 
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Figure 4: Real parts of Floquet multipliers for the fundamental breather at e = 0.01 near the 
bifurcation for the main branch (left) and side branches (right). Blue and red dots correspond to 
the Floquet multipliers with unitary and non-unitary absolute values respectively. 




T, Ji T,n 

Figure 5: Real parts of Floquet multipliers for the two-site breather with a hole at e = 0.01 near 
the bifurcation for the main branch (left) and side branches (right). Blue and red dots correspond 
to the Floquet multipliers with unitary and non-unitary absolute values respectively. 
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Figure 2] shows the real part of Floquet multipHers versus the breather's period for the fun- 
damental breather (left) and the new solution branches (right) bifurcating from the fundamental 
breather due to the pitchfork bifurcation. Because Floquet multipliers are on the unit circle for all 
periods below the bifurcation value, the fundamental breather remains stable for these periods, 
in agreement with Theorem [2l Once the the bifurcation occurs, one of the Floquet multiplier 
becomes real and unstable (outside the unit circle). Two new stable solutions appear during 
the bifurcation, which have the identical Floquet multipliers. These solutions become unstable 
for periods slightly larger than the bifurcation value, because of the period-doubling bifurcation 
associated with Floquet multipliers at —1. 

We perform similar computations for the two-site breather with the central hole (j44p . Figure [5] 
shows that at the coupling e = 0.01 the breather is unstable for periods 2tt <T < tJ*^^ and stable 
for periods T > tJ*^^ with tJ*^^ w 5.4257r. This can be compared using the change of stability 
predicted by Theorem [2l According to equation (f36|) , K2 changes sign from positive to negative 
at ri°^ 5.4767r. Since T'{E) is positive for the soft potential, Theorem [2] predicts that in the 
anti-continuum limit the two-site breather is unstable for T < T^^^ and stable for T > T^^\ This 
change of stability agrees with Figure [5] where we note that |Ti'^ - ri°^| ^ O.OSvr at e = 0.01. 

At T w S.Gvr and T ~ S.Tvr, two bifurcations occur and unstable multipliers bifurcate from 
the unit multiplier for larger values of T. The behavior of Floquet multipliers is similar to the 
one on Figure IHleft) and it marks two pitchfork bifurcations for the two-site breather with the 
hole. One bifurcation is visible on Figure [2] in the space of symmetric two-site breathers with 
n_i(t) = ui{t). The other bifurcation is observed in the space of asymmetric two-site breathers 
with n_i(t) 7^ ui{t). We display the two pitchfork bifurcations on the top panel of Figure [6j One 
can see for the second bifurcation that the value of oq is the same for both breathers splitting of 
the main solution branch. Although the values of a_i and ai look same for the second bifurcation, 
dashed and dotted lines indicate that ai is greater than a_i at one asymmetric branch and vice 
verse for the other one. The bottom panels of Figure [6] show the asymmetric breathers with period 
T = 5.757r. 

5.2 Five-site model 

We can now truncate the discrete Klein-Gordon equation ([T]) at five lattice sites, e.g. at n E 
{—2,-1,0,1,2}. The fundamental breather (j43p and the breather with a central hole are 
continued in the five-site lattice subject to the symmetry conditions Un{t) = U-n{t) for n = 1,2. 
We would like to illustrate that increasing the size of the lattice does not qualitatively change the 
previous existence and stability results, in particular, the properties of the pitchfork bifurcations. 

Figure [7] gives analogues of Figure [2] for the fundamental breather and the breather with a 
hole. The associated Floquet multipliers are shown on Figure [H We can see that both existence 
and stability results are analogous between the three-site and five-site lattices. 
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Figure 6: Top: The initial displacements a_i, gq, and ai for the T-periodic breather with a hole 
on the three-site lattice with e = 0.01. Bottom: Asymmetric breathers with period T = 5.757r on 
the three-site lattice with e = 0.01. 




a a a 

1 2 

Figure 7: Top: The initial displacements ao, ai, and a2 for the T-periodic fundamental breather 
of the five-site lattice with e = 0.01. Bottom: The same for the two-site breather with a hole. 
The dotted lines correspond to the T-periodic solutions to equation (jH]). 
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Figure 8: Top: Real parts of Floquet multipliers for the fundamental breather near the bifurcation 
for the main branch (left) and side branches (right). Bottom: The same for the two-site breather 
with a hole. Blue and red dots correspond to the Floquet multipliers with unitary and non-unitary 
absolute values respectively. Blue and red dots correspond to the Floquet multipliers with unitary 
and non-unitary absolute values respectively. 
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6 Symmetry-breaking bifurcation near 1:3 resonance 



We study here the symmetry-breaking bifurcation of the fundamental breather. This bifurcation 
ihustrated on Figure [3] occurs for soft potentials near the point of 1:3 resonance, when the period 
T is close to Gvr. For simplicity, we consider the soft quartic potential with V'{u) = u — . 

When u = 0^*^^ is the fundamental breather and T ^ 27rn is fixed, Theorem [1] and Lemma [1] 
imply that 

' UQ{t) = ip{t)-2eMt) + Oj,2^^(o,T)(e'), 
u±i{t) = +eMt) + Oj,2^^(o,T)(e'), (45) 

u±n{t) = + Oj^2^^(o,T)(e^), n>2, 

where ip can be expanded in the Fourier series, 

V9(t)= CniT) COS (46) 



dd 



and the Fourier coefficients {c,i(r)}„gN^j^ are uniquely determined by the period T. The correc- 
tion terms ifi and ipi are determined by the solution of the linear inhomogeneous equations (fT2] 
and (|13p . in particular, we have 



It is clear that if 27r < T < 67r is fixed, then ||ti±i||j?2j,j.(o,T) = ^(^) ^^'^ the cubic term is 
neglected at the order 0(e), where the linear inhomogeneous equation (jl2p arises. 

Near the resonant period T = Gvr, the norm ||it±i ||_f/2^^(o,T) is much larger than 0{e) and the 
cubic term Uj_i must be incorporated at the leading order of the asymptotic approximation. Let 
us introduce the scaling transformation, 

^= 1+^62/3 ' ^ = (1 + -^^'^')*' n„(t) = (l + <5e2/3)C/„(r), (48) 

where (5 is a new parameter, which is assumed to be e- independent. The discrete Klein-Gordon 
equation ([1]) can be rewritten in new variables as follows, 

Un + Un-U^ = fJ.Un + 7iUn+l + Un-l), H G Z, (49) 

where 

T-periodic solutions in variables {n„(t)}„gz become now 67r-periodic solutions in variables 
{UniT)}nez- To reduce the system of Klein-Gordon equation (j49p to a perturbed Duffing oscilla- 
tor, we consider the fundamental breather near 1:3 resonance and prove the following reduction. 
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Lemma 4 Let 5 be fixed independently of small e > 0. Let Uq be represented by the Fourier 
series, 

Uo{r)= Yl ^ncos(^), (51) 

such that the coefficients {6n}nGNodd converge to {c„(67r)}„gNodd as e ^ 0, where {c„(r)}„gNodd 
are uniquely defined by the Fourier series ^6\ )- Let a{6) be a root of the cubic equation, 

-26a{6) -ja^{6) =C3{6tt). (52) 

Assuming that 86 + 9a^((5) ^ and 6^0, there exists a solution of the discrete Klein-Gordon 
equation d^g] ) such that 

U±i{t) = ei/3a(5)cos(r) + e2/3_8|W_,os(^) + 0^,2^^.(0,6.) (e), 
U±2{r) = _e2/3f^eos(r) + OH|„(o,6.)(e), (53) 

U±n{T) = + C'//2^^(0,67r)(e), > 3. 

Proof. We shall consider the substitution, 

C/„(r) = ei/3a„cos(T) + K(t), n e N, (54) 
where a sequence {a„}„gN satisfies the system of difference equations, 

--a^-/ie"^/^a„-7e~2/3(a„+i+a„_i) =7e-^635n,i, n G N, (55) 
and the perturbation term {14i(''")}nGN satisfies the system of perturbed Duffing oscillators, 
Vn + Vn- - 3eV3a„ cos(t)V;2 - 3e2/3a2 cos2(t)K - - 7(K+i + K.-i) 

= iKi EfceNodd\{3} cos (^) + leal cos(3r), n E N. (56) 



Here we shall use oq = and Vq = in equations (j55p and (j56p for n = 1. Note that the 
resonance term cos(r) near the zero equilibrium has been removed from the series in (f56|) by the 
substitution (j54p . Recall that /i = ©(e^/^) and 7 = ©(e) as e — t- 0. Provided that the sequence 
{onlnGN is bounded as e — )• 0, the Implicit Function Theorem yields a unique continuation of a 
small Gvr-periodic solution with respect to parameter e > 0, which satisfies ||V||p(-p^ j|,2^^(o,T)) — 
for some C > 0. 

To solve the system of difference equations (|55l) for fixed (5 G M and small e > 0, we note that 
_^ 2(5 and je~^ — t- 1 as e — ?• 0. The assumption of the lemma states that 63 — )■ C3(67r) as 
e — )• 0. We shall consider a root of system (|55l) at e = given by a„ = a{5)5n,i, where a{5) is 
a root of the cubic equation (j52p . Under the conditions 85 + 9a'^{5) 7^ and 6^0, the Implicit 
Function Theorem yields a unique continuous of this root with respect to parameter e > 0, which 
satisfies 

a, = -eV3^ + 0(^2/3), 

an = + 0(e2/3), n>3. 

Substitution of this expansion to (|54l) yields the assertion of the lemma. ■ 
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Figure 9: Roots of equation (j52p . 

Remark 3 The condition 85 + 9a^{5) = implies bifurcations among the roots of the cubic 
equation \50fi . e.g., the fold bifurcation, when two roots coallesce and disappear across a bifurcation 
value of 5. The condition 5 = does not lead to new bifurcations but implies that the values of 
an for n > 2 are no longer as small as 0{€^^^). Refined scaling shows that if S = 0, then 
ai = a(0) + 0{e^/^), oa = 0{e^/^), and a„ = ©(e^/^?^^ ^ > 3^ ^j^^^^ ^^0) is a unique root of (5^ 
for 6 = 0. 

Roots of the cubic equation (I52p are shown on Figure [H One root continues across 6 = 
and the other two roots bifurcate for (5 < by means of a saddle-node bifurcation. Since (5 < 
for small e > correspond to T > 6ir, this saddle-node bifurcation gives a birth of two periodic 
solutions with the opposite sign in amplitude ai. This bifurcation is observed on Figure [2] (right), 
one of the two new solutions still satisfies the asymptotic representation (I45p as e — >• for fixed 
T > Qtt. 

We shall now focus on the root of the cubic equation (I52p that continues across 6 = 0. 
Substituting the asymptotic expansion (j53p into the discrete Klein-Gordon equation (j49p for 
n = 0, we obtain the perturbed normal form for 1:3 resonance, 

Uo + Uo-Ul = fiUo + i^cos(r) + 0^2^^(0,5.) (e'/'), (57) 

where u = 2-fe^/^a{6). Recah that if 5 e M is fixed, then = ©(e^/^) and u = 0{e^/^) as e ^ 0. 
We assume that a{6) 7^ 0, which is equivalent to the assumption that C3(67r) 7^ in the Fourier 
series (I46p . 

We are interested in Gvr-periodic solutions of the perturbed normal form ([57|) in the limit of 
small fi and u. Since the remainder terms are small as e — )• and the persistence analysis is rather 
straightforward, we obtain main results by studying the truncated normal form for 1:3 resonance, 

U + U -U^ = fiU + iycos{T). (58) 

A family of Gvr-periodic solutions with even symmetry exists in equation ([58|) for all small values 
of n and v. 




25 



Lemma 5 There are > 0, > 0, and C > such that for all ^ € (— MO) A'o) o-'iT'd v £ (— z^Oi ^^o); 
the normal form equation l[58\) admits a unique 6ir-periodic solution U^^u £ //g(0, Gvr) satisfying 

\\U^,.-v\\H^^^^<C{\fi\ + \u\). (59) 
Moreover, the map M x M 9 (/x, i^) h-^ [/^_^ G F^(0, 67r) is C°° (/x, i^) = (0, 0). 

Proof. The proof follows by the Lyapunov-Schmidt reduction arguments. For = and 
small n G (— //o,//o)j there exists a unique Gvr-periodic solution U^^ with even symmetry, which 
is 0{fi)-close to if in the F2^,(0,67r) norm. Because the Duffing oscillator is non-degenerate, the 
Jacobian operator L^^o has a one-dimensional kernel spanned by the odd function C/^,0) where 

L^,. = a,2 + l-/x-3[/2,(t). (60) 

By the symmetry arguments (f//i,07 cos(-))j;^2^^(-q g^) = 0, and the unique solution with even sym- 
metry persists for small v G {—i'q,i'q). ■ 



By the perturbation theory arguments, the kernel of the Jacobian operator ,^ is empty for 
small /i and v provided that v ^ 0. Indeed, expanding the solution of Lemma [5] in power series 
in fi and u, we obtain 

Uf,,u = ^ + fJ-L'^f + yL~^ cos(-) + OH2^^{ofi7T){lJ''^, 1^^), (61) 



where is the operator in (|13p . Although has a one-dimensional kernel spanned by (p, 
this eigenfunction is odd in r, whereas (p and cos(-) are defined in the space of even functions. 
Expanding potentials of the operator Lfj_u, we obtain 

Lt,,,yU^,u = z^sin(-) + C'//2^^(o,67r)(^^^^^^)• (62) 

We note that 

((^,COs(-))L2„(0,67r) = (<^,sin(-))L2^^(o,67r) / 

if C3(67r) 7^ 0, where C3{T) is defined by the Fourier series (|^6|) . By the perturbation theory, the 
kernel of L^^,y is empty for small i/ G (—1^0, z^o)- 

If the linearization operator L^^j^ becomes non-invertible along the curve v = z^*(^) of the 
codimension one bifurcation, the symmetry-breaking bifurcation occurs at = i^*(^). Figure [TOl 
(left) shows the behavior of Floquet multipliers of equation L^^^^W = with respect to parameter 
at /i = 0. We can see from this picture that the bifurcation occurs at ~ 0.00015. 

The right panel of this figure gives the dependence of the bifurcation value z/* from /x. Using 
the formula for ^ in (|50p . we obtain 

r = 6.^^. 

As the coupling constant e goes to zero, so does parameter u. As shown on Figure [TOl (right), 
parameter /i at the bifurcation curve goes to negative infinity as 1/ — )• 0. This means that the 
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Figure 10: Left: Floquet multipliers of equation Lfj_^^W = 0. Blue and red dots correspond to the 
Floquet multipliers with unitary and non- unitary absolute values respectively. Right: Parameter 
v versus /x at the symmetry-breaking bifurcation. 



¥ = 10"" ¥ = 2x10"" 




Figure 11: Solutions with period T = Gtt to equation (|58p at n = before (left) and after (right) 
the symmetry-breaking bifurcation. 

closer we get to the anti-continuum limit, the further away from Gvr moves the pitchfork bifurcation 
period T. 

Figure [TT] shows one solution of Lemma [5] for < < z^*(//) and three solutions for > i^*(//), 
where = 0. This behavior resembles the pitchfork bifurcation shown on Figure [3j 

Figure [T2] transfers the behavior of Figures [TO] and [TT] to parameters T, e, and = Uo{0). 
The dashed line on the left panel shows the dependence of period T at the pitchfork bifurcation 
versus e for the full system (I42p . We can see that the truncated normal form equation (I58p allows 
us to capture the main features of the symmetry-breaking bifurcations in the full system (j42p . 

7 Conclusion 

We have considered existence and stability of multi-site breathers in the Klein-Gordon lattices 
with harmonic lattice coupling. We have described explicitly how the stability or instability of a 
multi-site breather depends on the phase difference and distance between the excited oscillators. 
It is instructive to compare our results to those obtained by Yoshimura [T6] for the lattices 
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1.01 



Figure 12: Left: Period T plotted against e at the symmetry-breaking bifurcation of the funda- 
mental breather modelled by equation ([58|) (solid line) and equation (|^2|) (dashed line). Right: 
Bifurcation diagram for the initial displacement no(0) = ao and period T in variables (I48p com- 
puted from the Gvr-periodic solution to equation (j58p . 

with purely anharnionic coupling, 

Un + Un±v}^ = e{Un+l " ""n)^ - e{Un - Un-l)^ , (63) 

where A; > 3 is an odd integer. Table II summarizes the result of [16j for stable configurations of 
two-site breathers from the configuration 

u(°)(t) = aj(p{t)ej +akip{t)ek, 

where N = \j — k\ > 1. Note that the original results of [16j were obtained for finite lattices 
with open boundary conditions but can be extrapolated to infinite lattices, which preserve the 
symmetry of the multi-site breathers. 





N odd 


N even 


hard potential 
V'{u) = u + u^ 


in-phase 


anti-phase 


soft potential 
V'{u) = u — u^ 


anti-phase 
in-phase 


anti-phase 



Table II: Stable two-site breathers in soft and hard potentials from |16| . 

Table II is to be compared with Table I summarizing our results. Note that Table I actually 
covers M-site breathers with equal distance N between the excited sites, whereas Table II only 
gives the results in the case M = 2. We have identical results for hard potentials and different 
results for soft potentials. First, spectral stability of a two-site breather in the anharmonic poten- 
tials is independent of its period of oscillations and is solely determined by its initial configuration. 
This is different from the transition from stable anti-phase to stable in-phase breathers for even N 
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in soft potentials. Second, both anti-phase and in-phase two-site breathers with odd N are stable 
in the anharmonic lattice. The surprising stability of in-phase breathers is explained by additional 
symmetries in the anharmonic potentials. The symmetries trap the unstable Floquet multipliers 
at the point fi = 1. Once the symmetries are broken (e.g., for even A^), these multipliers split 
along the real axis and the in-phase two-site breather becomes unstable in soft potentials. 

We have also illustrated bifurcations of breathers near the point of 1 : 3 resonance. It is 
important to note that a similar behavior is observed near points of 1 : A: resonance, with k being 
an odd natural number. For the non-resonant resonant periods, a breather has large amplitudes 
on excited sites and small amplitudes on the other sites. As we increase the breather's period 
approaching a resonant point T = Ink for odd A;, the amplitudes at all sites become large, a cascade 
of pitchfork bifurcations occurs for these breathers, and families of these breathers deviate from 
the one prescribed by the anti-continuum limit. However, due to the saddle-node bifurcation, 
another family of breathers satisfying Theorem [T] emerges for periods just above the resonance 
value. The period-amplitude curves, similar to those on Figure [21 start to look like trees with 
branches at all resonant points T = Ink for odd k. In the anti-continuum limit, the gaps at 
the period-amplitude curves vanish while the points of the symmetry-breaking bifurcations go 
to infinity. The period-amplitude curves turn into those for the set of uncoupled anharmonic 
oscillators. 
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